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1.  Introduction 


This  is  one  of  a  series  of  reports  that  describes  the  Naval  Postgraduate 
School's  (NPS)  approach  to  marine  atmospheric  bo undry  layer  (MAHL)  modeling. 
The  basics  of  the  approach  and  the  status  of  modeling  and  parameterization 
of  the  pertinent  physical  processes  are  given  in  Fairall  et  al.  (1981).  The 
utility  of  the  model  for  tactical  use  and  initial  validation  of  the  model 
are  given  in  Davidson  et  al.  (1982).  In  this  report  we  describe  the  use  of 
the  model  for  analysis  and  prediction  of  the  aerosol  content  of  the  MABL. 

Interest  in  marine  aerosols  has  increased  recently  because  of  their  con¬ 
tribution  to  the  scattering  and  absorption  of  light.  Estimating  these  in¬ 
fluences  chi  electro-optical  (EO)  system  performance  has  been  emphasized  in 
several  studies  such  as  Bamhardt  and  Streete  (1970).  An  aerosol  size  of 
interest  is  that  associated  with  locally  generated  sea  salt  because  of  its 
effects  on  IR  as  well  as  visible  wavelengths.  Measured  distributions  in  the 
sea  salt  size  range  show  several  orders  of  magnitude  variability,  and  consi¬ 
derable  effort  has  been  expended  to  normalize  these  measurements  by  correct¬ 
ing  appropriately  for  relative  humidity  and  wind  speed.  These  two  meteoro¬ 
logical  quantities  are  considered  because  of  their  role  in  generation 
(wind),  transport  (wind),  and  growth  (relative  humidity)  of  aerosols. 

Both  empirical  and  theoretical  bases  exist  for  formulating  expressions 
for  equilibrium  aerosol  distributions.  Past  quantifications  of  the  depen¬ 
dence  of  equilibrium  distributions  on  relative  humidity  were  made  by  Fitz¬ 
gerald  (1975).  Dependence  on  turbulent  transport  was  considered  by  Tobc. 
(1965),  and  dependence  on  surface  generation  by  Chaen  (1973).  Ho  ;ver , 
aerosol  models  in  use  today  utilize  only  parameterizations  of  the  eff  - 
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of  relative  humidity  and  wind  speed  an  the  equilibrium  aerosol  distribu¬ 
tions  (e.g..  Wells  et  al,  1977).  Recent  evaluations  have  shown  that  these 
models  are  limited  to  mean  (that  is,  the  average  aerosol  density  encoun¬ 
tered  at  a  given  wind  speed  and  hunidity)  distributions  (Ehirall  et  al, 
1982).  These  models  are  inherently  limited  because  sane  processes  in  the 
atmospheric  mixed  layer  which  affect  aerosol  concentrations  are  not  con¬ 
sidered  . 

In  this  study  we  will  present  an  approach  for  including  meteorological 
descriptions  which  encompass  processes  and  structures  of  the  whole  marine 
atmospheric  boundary  layer .  The  top  of  the  boundary  layer  is  capped  by  the 
marine  inversion  where  entrainment  of  overlying  air  takes  place .  Because 
entrainment  mixes  clear  (non -marine)  air  into  the  marine  layer  this  process 
is  as  important  as  surface  layer  aerosol  fluxes  in  determining  the  equilib¬ 
rium  concentration.  The  primary  goal  of  our  examination  will  be  to  incor¬ 
porate  in  aerosol  descriptions  recently  established  features  of  the  inver¬ 
sion  capped  marine  boundary  layer.  The  features  wall  be  presented  in  terms 
of  an  integrated  slab  model .  Such  a  model  is  especially  desirable  since 
the  only  input  parameters  needed  are  routine  meteorological  observations . 


2.  Description  and  Evaluation  of  Existing  Aerosol  Models 

A  representative  example  of  current  models  for  estimating  equilibriim 
aerosol  distributions  is  tlat  formulated  by  Wells  et  al  (1977)  .  Their  for¬ 
mulation  vas  modified  by  investigators  in  the  tfevy  Electro-Optical/ Nfe teo- 
rology  Program  (Hughes,  1980).  The  number  density  spectrum,  n(r),  of  the 
sea  salt  component  is  described  in  the  modified  version  as 

-Z/h  F~8.5(r/a)Y 

n(r)  =  (r/a)  1.62  (C^-K^v6)/ Ffe  ,  (1) 

vhere  r  =  the  particle  radius  in  um 

u  =  the  wind  speed 
v  =  0.  5  for  u  <_  4  Tty's  , 
v  =  u  -  3  •  5  for  u  >  4  m/  s  , 

F  =  1  +  (v/60)3  , 

y  =  0.384  -  0.00293  v1-25  , 

Z  =  height  above  sea  surface,  m  , 
hQ  =  scale  height,  m  (800  m  for  Z  <  1000  m)  , 
a  =  0.81  exp  (0.066S/(1.058  -  S))  , 

S  =  H/100,  (H  is  relative  humidity  in  per  cent)  . 

The  other  constants  are 


v(ir\ /  s) 

C1 

C2 

6 

v  <_  7 

350 

1000 

1.15 

A 

> 

0 

6900 

0-29 

The  relative  humidity  growth  factor,  a,  has  the  form  sugqested  by  Fitz¬ 
gerald  (1975).  The  height  dependence  is  exponential,  with  the  scaling 
height,  h  ,  being  a  function  of  wind  speed  as  suggested  by  Tbba  (1901 
The  leading  term  (C^+C^  v  )  corresponds  to  local  generation  and  las  a  , 
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law  dependence  on  wind  speed.  This  particular  version  is  referred  to  as  the 
VMK  model. 

Examination  of  a  considerable  amount  of  data  obtained  in  the  Northeast 
Atlantic  and  the  Eastern  Pacific  indicates  that  Eqn  1  is  inccmplete.  Fig¬ 
ure  1  shows  the  height  dependence  of  total  aerosol  volume  frcm  a  sample  set 
of  Eastern  Pacific  data.  The  three  aerosol  profiles  correspond  to  (1)  the 
observed  sea  salt  aerosol  volume,  V,  (2)  the  observed  sea  salt  aerosol  vol¬ 
ume  adjusted  to  80%  relative  humidity  equilibrium  sizes,  V  ,  and  (3)  the 
VMK  predicted  volume  adjusted  to  80%  (open  circles).  These  data  were  ob¬ 
tained  under  forced  convective  conditions  and  with  active  local  production. 
In  order  to  examine  only  the  height  dependence  issue,  observed  and  model 
values  have  been  matched  at  the  surface  to  r amove  the  production  influence. 
It  is  clear  from  the  figure  that  within  the  mixed  layer  the  observed  de¬ 
crease  of  aerosol  volume  with  height  is  less  than  that  predicted  by  Bqn  1. 
The  surface  generated  aerosols  appear  to  be  well-mixed  belcw  the  inversion 
when  normalized  to  remove  the  influence  of  relative  humidity.  This  has  also 
been  observed  by  other  investigators  (Blanchard  and  Woodcock,  1980;  Hughes, 
1980;  Johnson  and  Hering,  1981).  Aerosol  volume  corrected  for  relative 
humidity  can  be  considered  to  be  reasonably  well  mixed  from  the  inversion 
down  to  approximately  the  10-meter  level  for  these  conditions. 

Evidence  supporting  our  assertion  that  existing  scaling  of  production 
effects  are  reasonable  only  for  mean  or  climatological  purposes  appears  in 
results  obtained  in  the  Northeast  Atlantic.  This  point  is  illustrated  in 
Fig.  2  where  normalized  measured  aerosol  volume  densities  at  r  =  5  um  ar  ; 
the  corresponding  prediction  frcm  Eqn  1  are  compared.  The  values  and  trends 
in  the  predicted  and  mean  results  are  in  reasonable  agreement.  However, 
the  standard  deviation  is  so  large  that,  at  a  given  time,  only  67%  of  the 
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Fig.  1  Height  dependence  of  aerosol  volune  and  relative  hunidity 

above  the  ocean.  The  curve  labeled  V  is  the  ambient  aero_»' 
3  3 

total  volune  ( urn  /cmJ),  the  curve  label  V  is  tiie  volune 

o 

corrected  to  00%  reference  humidity.  The  open  circles  re 
present  the  equivalent  VQ  fran  the  VMK  model  (Bq.  1) 


observed  aerosol  distributions  will  be  within  a  factor  of  three  of  the 
average.  The  point  is  that  no  matter  how  accurately  the  model  predicts  the 
average  aerosol  density  at  a  given  wind  speed  and  relative  humidity,  the 
factor  of  three  for  the  FMS  variation  cannot  be  eliminated  without  consi¬ 
dering  more  parameters  than  instantaneous  wind  speed  and  relative  hunidity. 
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3.  Boundary  Layer  Model 

The  prime  feature  characterizing  the  marine  atmospheric  boundary  layer 
is  that  it  is  ccanvectively  mixed,  and  hence  hcmogeneous,  up  to  a  height  h, 
where  it  is  capped  by  an  inversion .  The  mixed  layer  is  cool  and  moist  rela¬ 
tive  to  the  overlying  stable  layer  (shown  in  Fig.  3). 

The  fundamental  assumption  of  existing  models  of  the  boundary  layer  is 
that  veil-mixed  properties  r  ana  in  veil  mixed  when  undergoing  evolutions  in 
time.  If  the  value  of  any  veil-mixed  property,  X,  were  to  change  solely 
because  of  turbulent  fluxes  at  the  top  and  bottom  boundaries  and  uniform 
horizontal  advection,  and  if  X  ranains  veil  mixed  after  the  change,  then  it 
follows  that  the  vertical  flux  profile,  <W'X‘>,  must  be  linear.  For  this 
case  the  appropriate  balance  expression  for  the  local  change  of  X  is 
(Lilly,  1968) 

dX/dt  +  V  *VX  =  (<W'X'>  -  <W* X 1  >  • ) /h  ,  (2) 

n  O  1 

where  V^’VX  is  the  horizontal  advecticn ,  the  brackets  represent  a  suitable 
average  for  the  responsible  turbulent  scales,  and  subscripts  o  and  i  refer 
to  surface  and  inversion  values,  respectively. 

Difficulties  in  making  measurements  of  the  fluxes  at  the  inversion  and 
evaluating  the  horizontal  advecticn  make  Ecp  2  impractical  for  single  loca¬ 
tion  observations  without  further  simplification  and  scaling.  The  horizon¬ 
tal  advecticn  term  is  excluded  in  the  following  discussion  not  because  it 
is  necessarily  negligible  but  simply  to  limit  the  scope  of  this  paper.  Hue 
surface  flux  is  easily  calculated  using  the  bulk  aerodynamic  method,  ibis 
leaves  the  flux  at  the  inversion  to  be  parameterized  in  order  for  Eqn  2  to 
be  a  usable  balance  expression  for  single  station  assessments . 

The  relationship  between  the  equivalent  flux  property  at  the  inversior 
<W'X'>^,  the  rate  of  change  of  the  mixed  layer  depth  (dly^t)  >  the  mean 
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vertical  motion,  W,  the  entrainment  velocity,  Wg ,  and  the  "jimp"  in  X,  AX, 
at  the  inversion  is : 

dh/dt  =  W  +  We  ,  (3a) 

<W'X'>i  =  -  We  AX  .  (3b) 

The  utility  of  these  expressions  depends  on  the  capability  to  estimate  the 
quantity  Wg  from  bulk  parameters.  An  estimate  of  the  <W'X'>^  term  in  Eqn  3 
can  be  made  if  changes  of  the  inversion  height  and  subsidence  rates  are 
known  (W  =  dh/dt  -  W)  .  As  such,  it  would  not  be  useful  for  prediction 
since  cne  uses  Wg  to  predict  dh/dt,  but  it  could  be  applied  to  analytically 
explain  contributions  to  observed  changes  in  X. 

All  integrated  schemes  relate  the  availability  and  consumption  of  tur¬ 
bulent  kinetic  energy  at  the  inversion  to  surface,  cloud  regime,  and  in¬ 
version  parameters .  Turbulent  kinetic  energy  available  for  entrainment  is 
derived  from  surface  layer  and  cloud  region  buoyancy  fluxes  and  mechanical 
turbulence  associated  with  surface  lajer  and  inversion  lajer  wind  shears. 
The  effectiveness  of  the  available  energy  in  driving  entrainment  depends  on 
the  static  stability  of  the  inversion  lajer. 

The  simplest  formulation  for  Wg  is  for  the  cloud  free  zero-order  model 
(Fig.  3)  where  the  static  stability  is  scaled  by  the  jump  in  virtual  poten¬ 
tial  temperature.  The  available  kinetic  energy  is  assumed  to  be  a  fraction 
of  the  surface  buoyancy  flux.  The  expression  for  Wg  is  (Lilly,  1968) 

We  =  f  <W,9v'>o/A0v  ,  (4) 

where  suggested  values  of  f  lave  ranged  from  .1  to  .3. 

The  complexity  of  parameter! zaticns  of  this  type  increases  when  cloud 
region  buoyancy  is  included  and  vhen  more  realistic  structures  are  used  for 
the  inversion  zone.  Deardorff  (1978)  and  Stage  and  Busirger  (1981)  are 
suggested  for  more  complete  descriptions  of  the  scaling  procedures .  In  all 
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cases,  'however ,  a  useful  scaling  velocity  for  the  entrainment  rate  is  the 
convective  mixing  velocity, 

w*  =  (9Q0Zi/T>l/3  .  (5) 

where  =  <W‘  e‘v>0,  =  (1  +  a)h  is  the  inversion  height,  T  the  absolute 

temperature,  and  g  the  acceleration  due  to  gravity.  The  equivalent  expres¬ 
sion  to  Eq.  4  las  the  form 

We/W*  =  fW*2/(g2iA0v/T)  .  (6) 
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.  Aerosol  Considerations  in  a  Well-Mixed  Model 


Ihe  first  step  in  applying  the  mixed  lajer  model  to  aerosol  descrip¬ 
tions  is  the  establishment  of  the  well-mixed  parameters.  In  the  case  of 
temperature  (clear  sky)  it  is  the  virtual  potential  temperature;  in  the 
case  of  humidity  it  is  the  vater  vapor  mixing  ratio.  Although  a  given 
aerosol  particle  may  vary  in  size  due  to  changes  in  ambient  relative 
humidity,  the  salt  particle  spectrum  is  conserved  (because  of  the  lower 
number  concentrations  of  sea  salt  aerosols,  coalescence  can  generally  be 
neglected).  Thus,  the  well-mixed  property  for  aerosols  is  the  dry  size 
spectrum  mixing  ratio,  or  equivalently,  the  aerosol  spectrum  mixing  ratio 
at  sane  reference  relative  humidity.  Vfe  will  consider  the  aerosol  volume 
spectrum 

V(r)  =  4/3it  r3n(r)  ,  (7) 

where  n(r)  =  dN/dr  is  the  number  density  spectrum. 

Vfe  shall  define  V(r)  as  the  volume  of  aerosol  particles  per  an3  per 
radius  increment  at  STP  and  at  the  reference  saturation  ratio,  SQ  =  0.8. 

At  sane  height  in  the  mixed  lajer,  vhere  the  local  ambient  air  density  is  p 
and  the  saturation  is  S,  the  aerosol  volume  spectrum  is 

V'(rs(S))  =  V(rQ)  g2(S)p/po  ,  (8) 

where  Pq  is  the  density  of  air  at  the  surface,  rQ  and  rg  the  particle  radii 
at  the  reference  and  ambient  saturation 

rs  =  ro^S)  '  (9) 

and  g(S)  is  the  humidity  growth  factor 

g(S)  =  .81  exp(0* 066  S/(1.058-S))  .  (10) 

Note  that  p/ PQ  is  included  to  maintain  a  constant  mixing  ratio. 

The  next  step  in  applying  the  mixed  la^er  model  to  aerosols  is  to 
recognize  that  all  aerosols  do  not  originate  locally.  Both  above  and  bv 


the  inversion,  aerosols  are  advected  into  the  local  region.  In  terms  of 
aerosol  density,  entrairment  acts  as  an  aerosol  flux  (out  of  the  boundary 
layer)  because  the  aerosol  concentrations  above  and  below  the  inversion  are 
different.  In  the  mixed  layer  model  the  entrairment  acts  on  the  jimp 
across  the  inversion  and  the  flux  of  a  scalar  property  is 

<W,X'>i  =  -WgAX  =  -WjXp-X^)  ,  (11) 

vhere  is  the  quantity  just  above  the  inversion  and  X^  is  the  quantity 
just  below  the  inversion. 

With  aerosols  we  find  it  useful  to  consider  two  components  as  dis¬ 
tinguishable  from  each  other.  The  total  aerosol  volume  is  the  sum  of  the 
surface  generated  salt  component,  V  ,  and  the  background  (continental) 
component,  Vc  (Figure  4);  then 


V  =  V  +  V 

m  sn  an 

V  =  V 

P  cp 


Z  <  h 
Z  >  h 


(12a) 

(12b) 


where  we  have  assumed  =  0. 

sp 

At  this  point  Eqn  3  can  be  applied  to  each  component 

hfdV^/dt)  =  <W'VS’>0  -  (W^JV^  ,  (13a) 

h(  a^/dt)  .  we(vcp-vcm)  -%  ,  (13b) 

where  we  have  neglected  horizontal  advection,  assumed  negligible  local 
production  for  the  continental  component  (<w'vc,>0  ~  0)  and  included  the 
Stokes  gravitational  fallout  term  (WU,  1979).  Ihe  fallout  rates  above 
and  below  the  inversion  are  different  because  of  the  change  in  the  aeros 
spectra  caused  by  the  humidity  growth  factor.  The  Stokes  velocity  is 
calculated  from 


\  =  2g(  Pw-p)ro2g2(S)/(9Ep)  ,  U4, 
where  pw  is  the  density  of  the  droplet,  g  the  gravitational  accelerat  t 
and  E  the  kinematic  viscosity  of  air.  Since  the  rate  of  change  of  tl. 
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Altitude 


continental 


Aerosol  volume  Mixing  Ratio 
at  Reference  Relative  Humidity 


Pig-  4  Haight  dependence  of  continental,  Vc,  and  sea-salt,  V  , 

aerosols  in  the  mixed  la^er  model.  As  shovn,  the  aerosol 
volunes  are  represented  at  fixed  hunidity 


continental  component  ( dV^/dt)  due  to  mixed-layer  processes  is  often  anall 
compared  to  the  advective  effects,  no  attempt  will  be  made  to  investigate 
Ecp  13b.  Ihe  actual  process  of  determining  the  separate  components  of  a 
measured  aerosol  spectrum  presently  relies  on  the  feet  tlat  virtually  all 
particles  with  radii  less  than  0.1  p  are  of  "continental"  origin  (Ehirall 


5-  Efarine  Boundary  Layer  Experiment 

The  previously  described  aerosol  balance  expression  will  be  evaluated 
with  aerosol  and  atmospheric  mixed  layer  measurements  made  off  the  U.  S. 

West  Cbast .  They  were  made  during  the  Cooperative  Experiment  on  West  Cbast 
Oceanography  and  Meteorology  (CEWCCM-78)  which  was  conducted  frcm  25  April 
to  23  May  1978  by  several  U. S.  tfevy  sponsored  groups.  Measurements  were 
made  frcm  the  research  vessel  ACANIA  and  several  shoreline  radiosonde  sta¬ 
tions.  The  general  area  of  the  experiment  is  shown  in  Eigure  5.  Descrip¬ 
tions  of  synoptic  scale  conditions  influencing  the  Los  Angeles  basin 
throughout  the  entire  CEWCCM-78  period  have  been  presented  by  Rosenthal  et 
al  (1979)  .  Ebr  this  evaluation  we  chose  a  period,  20-21  ftey,  toward  the 
end  of  the  experiment. 

The  20-21  May  period  was  one  of  increasingly  maritime  conditions  in  the 
Los  Angeles  basin  due  to  onshore  flow  caused  by  the  developnent  of  a  ther¬ 
mal  trough  over  Southern  California .  This  thermal  trough  developed  in 
conjunction  wdth  the  intensification  of  a  Ehcific  high  located  west  of 
Washington  state . 

Mixed  la>er  changes  during  this  period  were  controlled  by  the  following 
factors:  1)  widespread  subsidence  was  occurring  over  a  uniform  mixed  layer, 
2)  prevailing  onshore  flow  reduced  the  effect  of  local  land/ sea  circula¬ 
tions,  3)  overcast  stratus-stratocunulus  contributed  to  the  entrainment  at 
the  top  of  the  la>er,  and  4)  advection  was  not  a  significant  factor  in 
local  changes  of  mixed  layer  depths . 

Deepening  of  the  marine,  or  mixed,  lajer  during  the  period  is  evid-  : 


in  an  acoustic  sounder  record  which  was  obtained  on  the  R/V  ACANIA  ana  i.rc, 
canposite  profiles  constructed  fran  radiosondes  la  inched  at  shoreline  and 
ship  locations  (Fig.  6).  The  acoustic  sounder  record  and  the  composite 
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Fig.  6a  Cfcmposite  radiosonde  measurements  of  0  at  5/19,  1700 
local  time 
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Fig.  6b  Cfcmposite  radiosonde  measurements  of  at  5/20,  0500 
local  time 
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Fig.  6d  Gtmposite  radiosonde  measurements  of  6v  at  5/21,  0500 
local  time  - 
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files  indicate  that  the  mixed  layer  depth  was  nearly  the  same  at  all  sta¬ 
tions  and  changed  uniformly. 


r 


Measurements  frcm  the  R/V  ACANIA  provided  information  on  the  surface 
layer  wind,  temperature,  humidity  and  aerosol  spectrum.  These  data  are 
summarized  in  Table  I  for  the  period  of  interest. 

Table  I 

Meteorological  surface  layer  data  for  the  CEWCCM-78  analysis 
period.  Tg  is  the  sea-surface  tenperature  and  T  the  air 
temperature.  These  are  four-hour  averages  of  half  hourly 
observations . 


Date 

Time 

u(  m/s) 

S 

TS(C) 

T(C) 

h(  m) 

5/20 

1300 

7.9 

.86 

14.0 

12.3 

425 

5/20 

1700 

8.7 

.82 

14.4 

12.6 

435 

5/21 

2100 

9.8 

.87 

13.7 

12.1 

385 

5/21 

0100 

9.2 

.88 

13.1 

12.0 

448 

5/21 

0500 

8.2 

.83 

13.0 

12.0 

615 

Practical  application  of  a  mixed  layer  model  requires  parameterization 
of  the  entrairment  rate,  Wg.  However,  in  Section  4  we  wish  to  test  only 
the  aerosol  part  of  the  model,  not  the  accuracy  of  various  entrairment  pre¬ 
dictions.  Therefore,  the  entrairment  rate  was  directly  estimated  frcm  the 
measured  evolution  of  the  wall-mixed  humidity  and  temperature  over  the  per¬ 
iod  using  the  NPS  boundary  layer  model.  An  entrairment  rate  was  selected 
that  reproduced  the  evolution  of  the  non-aerosol  mean  variables  (humidity 
and  temperature)  and  then  that  entrairment  rate  was  applied  to  the  aerosol 
evolution. 

The  vertical  motion,  W,  was  calculated  using  the  known  changes  in  h  and 
the  values  to  Wg  in  Eqn  3a.  The  divergence,  D,  was  estimated  as  D  =  -W/h 
(Table  II). 
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liable  II 


Ehtrairment  velocity,  We<  and  divergence  for  the  analysis 
period . 

Period  Wg ( an/ sec )  Divergence  (sec  1 ) 


5/20-0500  to 

0.32 

7.5  x  10-6 

5/20-1700 

5/20-1700  to 

0.32 

1.3  x  10"6 

5/20-2400 

5/21-0000  to 

0.38 

1.3  x  10-6 

5/21-0500 
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6.  Aerosol  Model  Test 


This  section  is  an  examination  of  the  ability  of  Efcjn  13a  to  predict 
evolutions  of  sea  salt  aerosol  spectra  given  ideal  estimates  of  the  meteo¬ 
rological  parameters,  consequently,  measured  (as  opposed  to  model  para¬ 
meterized  or  predicted)  values  are  used  for  non-aerosol  parameters  wherever 
possible.  In  particular,  the  wind  speed,  relative  humidity,  entrainment 
rate  and  boundary  layer  height  for  the  analysis  will  be  based  on  actual  ob¬ 
servations,  not  model  predictions, 
a.  Tterosol  Surface  Elux 

The  evolution  of  the  sea  salt  component  of  the  aerosols  is  de¬ 
scribed  by  Eqi  13a  where  h,  and  Wg  are  given .  The  only  other  unknown 
in  the  equation  is  the  surface  aerosol  flux  spectrum  <W'VS'>0.  This  quan¬ 
tity  has  been  determined  previously  from  data  collected  in  the  Northeast 
Atlantic  during  the  JASIN  experiment,  ( Ehirall  et  al,  1982)  and  is  given  in 
Thble  III  as  a  function  of  wind  speed . 

Table  III 

2  2 

Sea  salt  aerosol  surface  volume  flux  in  yin  /cm  / s,  as  a  func¬ 
tion  of  wind  speed  (m/s)  and  particle  radius  ( um)  at  S  =  0-8. 


u 

r  =  0.8 

2 

5 

10 

15 

6 

1.3 

1.1 

2.5 

1.0 

0-3 

9 

4.5 

3. 1 

4.2 

3.3 

2.3 

11 

8.2 

7.7 

11.0 

21.0 

27.0 

13 

9.1 

9.2 

17.0 

49.0 

48.0 

15 

11.0 

10.0 

19.0 

72.0 

140.0 

18 

17.0 

11.0 

24.0 

92.0 

180.0 

b.  Initialization 

Eqa  13a  is  a  rate  equation  that  relates  dV  /dt  to  V  .  Cbnse- 

sm'  am 

quently  the  model  is  used  to  predict  the  evolution  of  subsequent  to 
sane  initial  time  for  vhich  values  of  the  relevant  parameters  are  know  or 


assumed.  In  this  example,  we  will  use  the  values  measured  at  the  beginning 

of  the  analysis  period.  If  measured  aerosol  values  are  not  available,  the 

model  can  be  initialized  with  the  equilibrium  spectrum  obtained  fran  Eqn 

13a  with  dV/dt  =  0: 
sm 

Vsm  (Equilibrium)  =  <W'Vs'>0/(We  +  W,^)  .  (15) 

c.  Results 

The  model  was  initialized  for  1300  PDT  on  5/20/78  with  the  measured 
spectrum.  The  evolution  of  the  size  spectrum  was  calculated  with  half-hour 
time  steps  at  five  particle  sizes  (r  =  0.8,  2,  5,  10  and  15  urn)  using  the 
data  given  in  Tables  I,  II  and  III.  The  surface  flux  used  at  each  wind 
speed  was  a  log  interpolation  of  the  values  in  Table  3  weighted  by  the 
velocity  category.  The  results  of  the  model  calculation  for  three  radii 
are  shown  in  Figure  7  and  compared  with  the  measured  data  and  the  VMK  model 
at  four-hour  intervals.  Each  of  the  three  results  plotted  is  normalized  by 
its  cwn  particular  value  at  the  start  of  the  analysis  period.  This  was  done 
to  avoid  biasing  the  comparison  in  favor  of  the  mixed  layer  model  which, 
since  it  is  initialized  with  the  measured  data,  should  tend  to  automatically 
predicrt  the  correct  results  more  accurately  (especially  in  the  earlier 
periods).  The  initial  and  final  measured  spectra  are  shewn  in  Figure  8 
along  with  the  mixed  layer  and  model  predictions. 

d.  Discussfon 

The  mixed  layer  aerosol  parameterization  did  a  credible  job  of  p. 
dieting  the  aerosol  spectum  16  hours  after  initialization.  A  loq-averaq  Li 
yielded  an  average  of  the  predicted  to  the  measured  spectra  of  0.92  x  a 
1.54  (1.54  =  e°  where  o  =  0.43  is  the  standard  deviation  of  the  lev,  > 
ratio  and  <x>x  or  %e°  =  <x>e±a  =  e<4n  x>"a).  In  this  case,  the 
also  performed  well  with  a  log-average  ratio  of  prediction  to  measui .  ■ 


Ratio 


13  17  21  1  5  17  21  1  5  17  21  1  5 

Hour  of  Day 

Fig.  7  Ratio  of  ambient  aerosol  volime  at  the  hour  indicated  to 

the  volune  at  the  beginning  of  the  period,  V  (  r ,  t)  /V_(  r ,  o)  . 
The  open  circles  are  the  mixed  layer  dynamic  model,  the 
X' s  are  the  measurements  and  the  solid  circles  are  the 
WMK  model.  The  mixed  layer  model  and  measurements  were 
normalized  by  the  initial  value  of  the  data.  The  VMK  ratio 
was  obtained  by  normalizing  with  the  WMK  model  value  obtained 
from  the  initial  wind  speed  and  relative  hunidity 
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ratios  of  1.69  x  or  t  1.6.  This  is  uncharacteristically  good  performance 
for  the  VMK  type  of  model.  In  practical  application,  the  mixed  layer  model 
will  probably  not  perform  this  well  since  several  important  model  para¬ 
meters  (Wg,  h,  u  and  S)  were  not  model  predicted  but  obtained  frcm  the  mea¬ 
surements.  This  issue  will  be  examined  in  section  8.  Other  sources  of 
error  (such  as  advection)  will  also  have  an  impact.  The  most  important 
parameter,  the  surface  flux  of  aerosols  (Table  III),  was  based  primarily  on 
a  single  data  set  obtained  in  the  North  Atlantic.  Further  refinements  in 
the  surface  flux  can,  at  least  potentially,  greatly  improve  the  model  accu¬ 
racy.  The  surface  flux  data  also  need  to  be  extended  to  larger  particle 
sizes. 
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7.  Full  Dynamic  Model  Test 

Eqn  13a  of  Section  4  was  tested  for  a  20-hour  time  period  using  mea¬ 
sured  or  directly  inferred  values  of  the  relevant  model  parameters.  The 
aerosol  model  spectrum  was  initialized  at  the  beginning  of  the  analysis 
period  with  the  measured  spectrum  This  section  describes  a  simulation 
where  a  complete  dynamical,  cloud- topped  boundary  layer  model  of  the  Dear- 
dorf f  (1976)  type  is  used  to  provide  the  necessary  parameters  for  Eqn  13a. 
a.  Procedure 

The  full  mixed  layer  model  was  run  in  three  12-hour  blocks  begin¬ 
ning  at  1700  local  time  on  5/19/78  with  subsequent  re-initialization  of  all 
meteorological  variables  at  the  0000  and  1200  GMT  radiosondes  (1700  and 
0500  local  time).  The  aerosol  spectrum  was  initialized  at  1700  on  5/19 
using  Eqn  15  and  was  not  subsequently  reset.  The  evolution  of  the  aerosol 
spectun  was  calculated  using  Eqn  13a  with  all  parameters  taken  frcm  the 
dynamic  model  output.  The  measured  aerosol  data  were  available  for  the 
last  20  hours  of  the  38-hour  model  simulation  (the  final  time  period  was 
run  for  14  hours  to  cover  the  last  two  hours  of  aerosol  data) .  Other 
details  of  relevance  are: 

1)  The  surface  fluxes  were  obtained  using  standard  bulk  aero¬ 
dynamic  parameterizations. 

2)  The  sea-surface  temperature  was  assigned  as  the  measured  valtv 
at  the  beginning  of  a  12-hour  forecast  period  and  'ssunod  to  \ 
constant  throughout  that  period. 

3)  The  measured  wind  speed  for  a  12-hour  period  was  simp]ifi>' 
two  or  three  continuous  linear  segments  in  an  effort  co 
late  an  accurate  wind  speed  forecast. 
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4)  The  subsidence  rate  for  the  start  of  the  test  vas  based  on  an 
assumed  climatological  divergence  of  D  =  1  x  10-^  sec-'*'.  The 
difference  between  the  predicted  and  measured  values  of  h  at 

-a  the  end  of  each  12-hour  period  was  used  to  re-estimate  the  ap- 

a  propria te  subsidence  rate  for  the  period,  this  new  subsidence 

value  was  then  applied  to  the  next  12- hour  period. 

5)  Mvection  of  temperature,  water  vapor,  mixed-layer  depth  and 
aerosols  were  assumed  to  be  negligible. 

b.  Aerosol  Results 

Mixed  layer  model  and  VMK  model  aerosol  spectra  were  calculated  as 
described  in  Sections  5.1  and  4.3.  The  results  are  shown  in  Fig.  9  for  the 
entire  38  hour  period.  Note  that  in  this  case  the  actual  spectral  values 
are  plotted  rather  than  the  model/ initial  ratios  shown  in  Fig.  7.  The 
mixed  layer  model  was  not  a  significant  improvement  over  the  VvMK  for  this 
data  set.  A  log-averaging  yielded  an  average  ratio  of  mixed  layer  model  to 
measured  aerosol  spectra  of  0.76  x  or  i  2.0  while  for  the  WMK  model  the 
ratio  was  0.87  x  or  j  2.2.  In  other  words,  the  mixed  layer  model  predic¬ 
tions  were  typically  within  a  factor  of  2. 0  while  the  VvMK  predictions  were 
typically  within  a  factor  of  2.2.  These  can  be  contrasted  with  the  results 
of  Section  4-4  where,  using  measured  parameters  and  spectral  initializa¬ 
tion,  the  mixed  layer  model  was  typically  within  a  factor  of  1.5  and  the 
VvMK  was  within  a  factor  of  1.6.  Thus,  a  more  realistic  application  of  the 
model  doubled  the  uncertainty  (note  that  a  factor  of  1.0  is  no  uncertain¬ 
ty)  . 

c.  Mixed  layer  Dynamics 

The  practical  application  of  the  full  dynamic  model  (Section  5)  led 
to  an  increase  of  uncertainty  in  the  aerosol  prediction  from  a  factor  of 
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Fig.  9  Temporal  evolution  of  the  aerosol  volime  at  r  =  0.8  um, 

5  um  and  15  um  ever  the  38  hour  period  beginning  at  1700 
local  time  on  5/19.  The  open  circles  are  the  full  mixed 
la>er  model,  the  solid  circles  are  the  WMK  model  and  the 
X*  s  are  the  data  • 
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1.5  to  a  factor  of  2.  The  primary  sources  of  the  increase  in  error  are: 

1)  Use  of  Eqn  15  for  initialization. 

2)  Madel  calculations  of  Wg. 

3)  Linear  smoothing  of  the  wind  speed  temporal  evolution. 

4)  Model  calculation  of  h. 

A  comparison  of  model  and  measured  values  of  Wgl  h  and  W  is  shown  in  Pig. 
10.  Recall  that  the  model  was  initialized  at  12  hour  intervals  based  on 
the  radiosondes  rather  than  surface  based  measurements  of  temperature  (e), 
hunidity  (q)  and  mixed  layer  depth  (h)  .  Whether  the  substantial  differences 
between  the  radiosonde  and  surface  based  data  are  due  to  measurement  inac¬ 
curacy  or  a  difference  in  sampling  is  not  known. 

The  calculation  of  evolutions  of  the  well  mixed  e  and  q  variables 
is  shown  in  Fig.  11.  Since  the  Deardorff  model  substantially  over-esti¬ 
mated  Wg  through  the  analysis  period,  the  model  predicted  too  much  warming. 
This  was  not  oily  because  of  the  entrainment  of  too  much  warm  air  but  also 
because  of  the  calculated  reduction  of  long  wave/ cloud- top  radiative  cool¬ 
ing  due  bo  overthinning  of  the  stratus  cloud  layer  by  the  model .  Another 
factor  in  the  prediction  of  9  and  q  was  the  overestimation  of  the  surface 
fluxes  caused  by  the  use  of  a  sea  surface  temperature  higher  than  the 
average  for  the  12  hour  period  (the  sea  surface  temperature  declined  by 
roughly  0.8  °C  during  each  12  hour  period) . 

The  tendency  of  the  Eeardorff  model  to  overestimate  Wg  for  thin  clouls 
has  been  discussed  by  Stage  and  Businger  (1981)  .  They  use  a  slightly  dif¬ 
ferent  parameterization  and  find  substantially  lower  values  of  Wg  if  Zc/h 
is  close  to  one  (Zc  is  the  lifting  condensation  level) .  A  comparison  of  W£ 
estimates  of  Deardorff  (1976)  and  Stage  and  Businger  (1981)  is  shown  in 


5/19  5/20  5/21 

Fig-  10  Ehtrainment  rate  ( W^)  and  mean  vertical  motion  (W)  are  on 
the  upper  scale  and  the  mixed  layer  depth  (h)  is  an  the 
lover  scale.  The  symbols  are  defined  in  the  graphs.  The 
discontinuities  in  the  model  results  are  due  to  re- initial 
zation  at  12  hour  intervals 
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fig.  11  Tanporal  evolution  of  the  wall  mixed  meteorological 

variables  6  and  q.  The  surface  data  were  measured  an  the 
R/V  ACANIA 
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T&ble  I V 


A  comparison  of  Deardorff  (1976),  Wg,  Stage  and  Businger  (1981), 
W  '  ,  entraiment  rates  for  the  initial  conditions  at  each  of  the 
12-hour  analysis  periods. 


Cate 

Time 

h(m) 

Zc(m) 

Zc/h 

W  (cm/s) 

W  '  (an/s) 

5/19 

1700 

310 

cio 

0.03 

0.86 

0.86 

5/20 

0500 

380 

270 

0.71 

0.84 

0.42 

5/20 

1700 

520 

340 

0.65 

1.43 

0.29 

Che  final  point  to  consider:  the  evolution  shown  in  Fig.  9  is  cal¬ 
culated  at  S  =  Sq  =  0.8.  Thus,  applications  requiring  the  ambient  spectra 
including  the  humidity  influence  (for  example,  optical  extinction)  are  sub¬ 
ject  to  additional  uncertainty  due  to  erroneous  model  predictions  of  S 
(Fig.  12)  and  therefore  g(S).  This  is  particularly  critical  if  S  >  0.9. 
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°  Model 
®  Radiosondes 
—  Surface  data 


Time(hrs) 
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Fig.  12  Similar  to  Figure  11  but  for  the  ambient  saturation 
ratio  S  at  Z  =  10  in 
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8.  Conclusions 


The  application  of  modem  mixed  layer  physics  to  the  structure  and  evo¬ 
lutions  of  atmospheric  aerosols  in  the  marine  regime  represents  a  philoso¬ 
phical  quantun  junp  over  the  primarily  empirical  formulations  of  the  VMK 
type.  This  was  accomplished  by  partitioning  the  aerosol  spectrum  into  con¬ 
tinental  and  locally  generated  components  and  by  writing  the  equations  in 
terms  of  aerosol  spectral  densities  transformed  to  a  predetermined,  fixed 
reference  relative  humidity.  Besides  the  wind  speed  and  relative  hunidity, 
the  new  formulation  requires  knowledge  or  specification  of  the  entrairment 
rate,  the  mixed  la^er  height  and  the  sea-surface  production  rate  of  drop¬ 
lets.  If  dynamic  effects  are  ignored,  the  model  can  predict  an  equilibrium 
aerosol  spectrum  (Eqn  15)  roughly  equivalent  to  that  of  the  WMK  model. 

The  aerosol  dynamic  model  was  tested  against  a  marine  data  set  in  two 
ways .  In  Section  6 ,  EJgn  13a  was  evaluated  under  optimum  computational  cir¬ 
cumstances  using  measurements  of  the  relevant  variables  wherever  possible. 
In  Section  7  ,  a  more  realistic  test  was  performed  where  a  full  mixed  la>er 
dynamic  model  evolution  was  used  to  simulate  an  actual  field  application  of 
the  aerosol  model .  In  both  cases  the  aerosol  mixed  layer  model  performed 
reasonably  well  but  it  did  not  convincingly  outperform  the  much  simpler  VMK 
model .  The  errors  in  prediction  of  the  aerosol  spectrum  and  the  ambient 
hunidity  represent  about  a  f a c tor-o f- three  uncertainty  in  the  aerosol  opti¬ 
cal  extinction  coefficient. 

Although  the  mixed  layer  model  was  only  a  modest  improvement  over 
VMK  model,  it  holds  considerable  potential  for  improvement.  The  inclur 
of  advective  effects  and  improvanents  in  the  entrainment  rate  and  at 
surface  production  flux  should  be  pursued.  Another  factor  to  oonsidf- 
the  nature  of  the  aerosol  production  used  here .  The  flux  data  given  . 


T&ble  III  are  actually  the  net  surface  flux  values  which  represent  the 
droplet  production  minus  losses  die  to  turbulent  deposition  under  average 
conditions .  A  modification  of  Ihble  III  to  represent  the  true  production 
and  the  inclusion  of  the  turbulent  deposition  velocity  (SI  inn  and  SI  inn, 
1980)  explicitly  in  Ecp  13,  is  presently  under  investigation. 
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